Critical behaviour of random fibers with mixed Weibull distribution 
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A random fiber bundle model with a mixed Weibull distribution is studied under the Global Load 
Sharing (GLS) scheme. The mixed model consists of two sets of fibers. The threshold strength of 
one set of fibers are randomly chosen from a Weibull distribution with a particular Weibull index, 
and another set of fibers with a different index. The mixing tunes the critical stress of the bundle 
and the variation of critical stress with the amount of mixing is determined using a probabilistic 
method where the external load is increased quasistatically. In a special case which we illustrate, 
the critical stress is found to vary linearly with the mixing parameter. The critical exponents and 
power law behaviour of burst avalanche size distribution is found to remain unaltered due to mixing. 



PACS numbers: 46.50. +a, 62.20.Mk, 64.60.Ht, 81.05.Ni 



I. INTRODUCTION 

Sudden catastrophic failure of structures due to unex- 
pected fracture of component materials is a concern and 
a challenging problem of physics as well as engineering. 
The dynamics of the failure of materials show interesting 
properties and hence there has been an enormous amount 
of study on the breakdown phenomena till nowi*^. The 
complexity involved in fracture processes can be suit- 
ably modelled by grossly simplified models. The simplest 
available model is Fiber Bundle Model(FBM)£»Mil. 

A FBM consists of N parallel fibers. The disorder of a 
real system is introduced in the fiber bundle in the form 
of random distribution of strength of each fiber taken 
from a probability density p(cr) and hence called Ran- 
dom Fiber Bundle Model (RFBM). The strength of each 
fiber is called its threshold strength. As a force F is ap- 
plied externally on a bundle of N fibers, a stress a = F/N 
develops on each of them. The fibers which have their 
threshold strength smaller than the stress generated, will 
break immediately. The next question that arises is the 
affect of breaking of these fibers on the remaining in- 
tact fibers i.e., one has to decide a load sharing rule. 
The two extreme cases of load sharing mechanisms are 
Global Load Sharing (GLS) and Local Load Sharing 
(LLS)±k±i^. In GLS, the stress of the broken fiber is 
equally distributed to the remaining intact fibers. This 
rule neglects local fluctuations in stress and therefore is 
effectively a mean field model with long range interac- 
tions among the elements of the system^. On the other 
hand, in LLS the stress of the broken fiber is given only 
to its nearest surviving neighbours. It is obvious that 
the actual breaking process involves a sharing rule which 
is in between GLS and LLS. Several studies have been 
made which considers a rule interpolating between GLS 
and LLSii^. 

For a given force F, some fibers break and they dis- 
tribute their load to the surviving fibers following a load 
sharing rule causing further failure and redistribution of 
stress. This process continues until all the remaining 
fibers have their threshold strength greater than the re- 
distributed stress acting on them. This corresponds to 



the fixed point of the dynamics of the system. As the 
applied force is increased on the system, more and more 
fibers break. An avalanche of size S is defined as the 
number of failed fibers between two successive external 
loadings. There exist a critical load (or stress a c ) be- 
yond which if the load is applied, complete failure of the 
system takes place. Most of the studies on FBM involve 
the determination of the critical stress a c and the inves- 
tigation of the type of phase transition from a state of 
partial failure to a state of complete failure. It has been 
shown that a bundle following GLS has a finite value of 
critical stress and belongs to a universality class with a 
specific set of critical exponent o 8 ' 16 whereas there is no 
finite critical stress a c at thermodynamic limit in the case 
of LLS in one dimensio n 10 : 18 . On the other hand, LLS on 
complex network has been shown to belong to the same 
universality class as that of GLS with the same critical 
exponents^. 

In this paper, we study a FBM with mixed fibers. 
Fibers have their threshold strength randomly chosen 
from Weibull distributions with two different index pa- 
rameters. The motivation here is to study the dynamics 
of random fibers in the presence of disorder caused due to 
mixing of two types of fibers with overlapping distribu- 
tion of threshold strengths. Moreover, the probabilistic 
method implemented to estimate the critical stress can 
be very easily used to study any type of mixed fiber bun- 
dles thus enabling one to put maximum disorder in the 
system being studied. 

Section II consists of description of the model high- 
lighting the method used. In section III we present the 
results. Section IV includes discussions and conclusions. 



II. THE MODEL 

In the Weibull distribution (WD) of threshold strength 
of fibers, the probability of failure of each element when 
a stress a is generated has a form 



P(a) 
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where cto is a reference strength and p is called the 
Weibull index. We consider a mixed RFBM with WD of 
threshold strength of fibers where strengths of fibers are 
randomly chosen from two different distributions charac- 
terised by different Weibull indices and study the critical 
behaviour of the model. A fraction x (henceforth called 
the mixing parameter) of fibers belong to the class A 
(WD with p — pi) and the remaining (1 — x) fraction be- 
long to the class B (WD with p — p 2 ) with the reference 
strength cr set equal to one in both the distributions. 

The probabilistic method introduced by Moreno, 
Gomez and Pacheco^ is extended to explore the critical 
behaviour of the above model. For the conventional WD 
(x = 1 or 0), let us consider a situation where the stress 
on each fiber increases from o\ to a 2 . The probability 
that a fiber randomly chosen from the WD survives from 
the load a\ but fails when the load is a 2l is given by 



P{<7 2 ) -Pfa) 
1-P(<7l) 
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Thus the probability that the chosen fiber that has sur- 
vived the load o\ also survives the (higher) load a 2 is 
<z(CTi,f2) = e~( <T2_<Tl ). The key point is that the force F 
on the bundle is increased quasistatically so that only the 
weakest fiber amongst the remaining intact fibers break. 
Thus, one needs to identify the weakest fiber, break that 
fiber and then calculate the number of remaining unbro- 
ken fibers using q(a±, <r 2 ) recursively updating the value 
of a 2 and o\ due to the load given away by the broken 
fibers, until no more failure occurs. As mentioned previ- 
ously, the dynamics of breaking will continue till the sys- 
tem reaches a fixed point. The process of slow increase 
of external load is carried on up to the critical stress a c . 
The method avoids the random averaging involved in a 
Monte Carlo simulation and hence turns out to be very 
useful in dealing with fluctuations. 

However, to deal with the 'mixed' random fiber bundle 
model, the above method needs to be generalised in an 
appropriate manner. The essential point is the fact that 
we need to keep track of number of unbroken fibers in 
each distribution separately. Let us assume that after 
a loading is done and a fixed point is reached, N kl and 
Nk 2 are the number of unbroken fibers corresponding to 
p = Pi and p — P2 distribution respectively, and a k is 
the stress per fiber at that instant. 

We define q\ and q 2 as follows: 

qi(<Ti,a 2 ) = e-^"-<^ 
q 2 ((Ti,a 2 ) = e-^ ~^ \ 

One needs to calculate the load (N^ + Nk 2 )ai that has 
to be applied to break one fiber. Here, 07 = min(<T; 1 , cr; 2 ) 
where 07 x (or ai 2 ) is the next weakest fiber in the p = 
px (or p = pi) distribution of strength of fibers and is 
obtained by the solution of the following equations: 



and N k2 
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FIG. 1: Variation of fraction of broken fibers with external 
load per fiber a for a RFBM with two distributions corre- 
sponding to pi = 2, p2 — 3, x = 0.5 and N=50000. 



which gives 
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(2) 
(3) 



The breaking of one fiber and the redistribution of its 
stress to all the remaining intact fibers causes some more 
failures. Let us assume that during this avalanche, at 
some point before the fixed point is reached, there are 
Nk-L and N k2 number of unbroken fibers belonging to the 
two distributions where each fiber is under a stress Ok. 
This stress causes some more failures and as a result N, 



fei 



and N u fibers are unbroken 

K2 



Let the new stress devel- 
oped be a k . The number of fibers which survive Ok and 
a k are obtained using the relation 

K= N 'kM^,a k ) (4) 
Nk 2 =N' k2 q 2 (^ k ,a k ) (5) 

The stress on each fiber is now equal to (N k2 + 

N k Ja' k /(N' k \ + N k J. Eq (4) and Eq. (5) are used again 
and again until a fixed point is reached. The fixed point 
condition is given by (N k +N k2 ) — (N k +N k ) < e where 
e is a small number (0.001). It should be noted that the 
critical behaviour does not depend on the choice of e. 
After the fixed point is reached, stress 07 is calculated 
once again as mentioned before and the whole process is 
repeated till complete failure occurs at a c . 

It should also be mentioned here that the critical stress 
<7 C can also be derived using directly the cumulative dis- 
tribution function for the mixed model as given in Eq. (6) 
in p(ai,<7 2 ). However in this case, the expression of 07, 
as defined above, turns out to be very complicated and 
it is difficult to arrive at a simple closed form as shown 
in Eq. (2) or (3). 
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III. RESULTS 

We present here the main results of a particular case, 
x = 0.5, p\ = 2 and pi — 3. Fig. 1 shows the fraction 
of total number of broken fibers as a function of applied 
stress a . The graph clearly shows the existence of a 
critical stress a c ~ 0.46 at which fraction of failed fibers 
increases rapidly and the bundle breaks down completely. 
The critical stress of the mixed bundle lies between that 
of the two pure bundles (for p = 2, a c = 0.42 and for p = 
3, <j c — 0.49 obtained using a c = ( / oe)~ 1 / p )£. Thus, the 
resulting critical stress of the mixed fiber bundle model 
can be tuned by varying the mixing parameter x. 

The mean avalanche size S of failure is defined as the 
total number of broken fibers between two successive 
loadings. It diverges near the critical point as (a c — er)~ 7 
with an exponent 7 = 1/2. Scaling behaviour of S is 
shown in Fig. 2. 

3.0 1 1 1 1 1 1 1 1 1 r 



2.5- 




-4 -3 
log 10 (aoc) 

FIG. 2: Scaling behaviour of mean avalanche size S as the 
critical point is reached for the mixed model. Also shown is a 
straight line (dotted) with slope (-1/2). Here, N=50000 and 
x = 0.5. 




FIG. 3: Avalanche size distribution for x=0.5, pi — 2, p2 = 3 
and N=50000. A straight line with slope -5/2 is also shown. 

The important feature associated with the failure pro- 
cess in RFBM is the power law behaviour of burst 



avalanche distribution of fibers. The probabilistic 
approach 9 implemented to determine a c turns out to 
be inappropriate in exploring the behaviour of avalanche 
size distribution. This is due to the fact that the mean 
avalanche sizes (S) obtained for different a are of frac- 
tional values which leads to a difficulty in calculating the 
distribution of avalanche sizes. We therefore use the stan- 
dard Monte-Carlo method along with the weakest fiber 
approach^ 7 , where the external load is increased by an 
amount sufficient to break the weakest intact fiber. The 
corresponding integer value of the size of an avalanche 
is denoted by A. For GLS, the distribution -D(A) of 
an avalanche of size A follows a power law D(A) oc A~^ 7 
where £ = 5/2 in the asymptotic limit 17 . Fig. 3 shows the 
avalanche size distribution for the present mixed RFBM 
obtained numerically with x = 0.5; a power law be- 
haviour with the same exponent 5/2 is clearly observed, 
confirming the mean field nature of the model. That the 
5/2 behaviour is expected even for a mixed RFBM for 
any x can be justified using the saddle point method ap- 
plicable in the limit of large Ail. In the present case 
(pi = 2 and p 2 = 3), the probability that a fiber will 
break when subjected to a stress a is 



P{a) = x[l - exp(-cr 2 )] + (1 - x)[l - cxp(-o- 3 
so that the density distribution becomes 
p(a) = 2xcrexp(-cr 2 ) + (1 - a;)3a 2 exp(-CT 3 ). 
The avalanche size distribution takes the form 



(6) 



(7) 
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l-P(a) 



exp( 
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l-P(a) 
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(8) 



where a* is the redistributed stress at the critical point 
at which the average applied force [F — 7Ver(l — P(a))] 
maximises, and p(a) and P{a) are as defined above. The 
function inside the square bracket has a maximum when 



l-P(a) 



1. 



(9) 



It should be noted that the above condition is satisfied 
when a — a* . Since the threshold distribution of fibers, 
(Eq. (7)), does not have any discontinuity, the saddle 
point integration of Eq. (8) (retaining the first order 
term in the expansion of the prefactor (1 — P(a) —<jp{<j)) 
around a — a*) yields the asymptotic behaviour -D(A) oc 
A- 5 / 2 . 

Let us now comment on the imminent failure 2 ^ be- 
haviour, when a fraction of weak fibers are already re- 
moved, and the distribution is close to the critical distri- 
bution i.e., the strength of the weakest intact fiber (<to) 
is close to the redistributed stress a* at the critical point. 
We shall consider the case when pi = 2 and pi = 3. The 
variation of D(A) with A is shown in Fig. 4 and as in the 
pure Weibull case, we observe a crossover from A -3 / 2 to 
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FIG. 4: The variation of D(A) with A when close to the crit- 
ical distribution for a mixed Weibull distribution with pi = 2, 
p2 — 3 and x = 0.5. A clear crossover from A _3//2 to A~ 5 ^ 2 
is observed with ao = 0.62. The slope of the thick line is -3/2 
and dotted line is -5/2. The inset shows the avalanche size 
distribution for the critical case along with a line of slope of 
-3/2. 

A~ 5 / 2 as A increases. For the critical distribution, how- 
ever, we observe a A -3 / 2 behaviour for the whole range 
of A (Inset Fig. 4). 

Although the power law behaviour of the avalanche 
size distribution near the critical distribution remain un- 
altered by mixing, one may ask the question about the 
behaviour of A c with x (where A c denotes the avalanche 
size at which crossover from 3/2 to 5/2 is observed). For 
the case of interest, p\ — 2 and p2 = 3, A c does not 
change appreciably with x. This is due to the fact that 
a* is almost constant as x is varied (explained later) . On 
the other hand, if we consider the case p\ = 1 and p2 — 3 
(Fig. 6), A c should decrease with x keeping er constant 
because as x is increased a* also increases taking the sys- 
tem away from critical distribution. 

The critical exponent 7 and the power-law behaviour 
of avalanche size distribution of the mixed model remain 
unchanged for any value of mixing parameter x. This 
supports the fact that the critical behaviour of a fiber 
bundle model is determined entirely by the load-sharing 
rule. 

However, the variation of <r c with the mixing parame- 
ter x shows a very interesting linear behaviour (Fig. 5). 
To justify the linear behaviour, we recall the probabil- 
ity distribution (Eq. (6)). A relation between the critical 
stress and a* can be obtained by calculating the applied 
force at a* where a* is a solution of Eq. 9. 

°c = ^ = cr* 2 p(a*). (10) 

For the conventional WD, a* = (-)?, which imme- 
diately reveals the fact that the difference in a* for the 
cases p = 2 and p — 3 is very small as compared to the 
change in corresponding a c values. To a good approxi- 
mation, one can set a* — constant = c, so that Eq.(10) 
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FIG. 5: Variation of critical stress a c with x. Slope of the 
straight line is 0.068. N=50000. 
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FIG. 6: Nonlinear variation of critical stress with x for p\ = \ 
and p2 — 3. 

combined with Eq. (7) yields 

a c = c 2 [2xcexp(— c 2 )] + (1 — x)3c 2 exp(— c 3 ) 

= xc 2 [2cexp(-c 2 ) - 3c 2 exp(-c 3 )] + 3c 4 exp(-c 3 ), 

which justifies the dominant linear dependence of a c on 
x for this particular case. It should also be emphasised 
that this linear relationship is a characteristic of situa- 
tions where the variation of a* is negligible (to the lowest 
order) as x is tuned from zero to one (e.g., the present 
case) and in general there is a non-linear relationship as 
shown in Fig. 6. 

In a recent worhSi, the shear failure of a glued inter- 
face has been studied using a simple beam model where 
beams (fibers) connect the two surfaces. The stretching 
and bending threshold strength of the beams, denoted 
by €1 and €2 respectively, are randomly distributed vari- 
ables satisfying a joint probability distribution function 
p(eij£2)- The mean field critical exponents arc obtained 
when the threshold distributions for bending and stretch- 
ing modes are independent and chosen from two differ- 
ent WD. In our model, we study fibers with threshold 
strength chosen from two Weibull distribution with dif- 
ferent Weibull indices and the critical exponents stick to 
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the mean field values also in our case. 

IV. CONCLUSION 

The critical stress of a mixed RFBM with Wcibull 
distribution and GLS is studied using a probabilistic 
approach where the external force is increased qua- 
sistatically at every step of loading^. The advantage 
of this method is that it does not require the process 
of random averaging which takes comparatively longer 
computational time. We obtain the variation of crit- 
ical stress with parameter x. This functional depen- 
dence of the critical stress on the characteristic quantity 
'the mixing parameter' is an important objective of this 
work. 

The critical behaviour of the mixed model namely the 
critical exponents and the power law behaviour of the 
burst avalanche distribution are same as mean field. In a 



mixed Weibull distribution, the threshold distribution of 
the two types of fibers are overlapping or the distribution 
is continuous. The presence of discontinuity modifies the 
avalanche size distribution for smaller Hence one 

expects the mean field (GLS) behaviour of the mixed 
model. The behaviour of the imminent failure shows a 
crossover in the avalanche size exponent from 5/2 to 3/2 
as the critical distribution is approached. In some distri- 
butions (depending upon p\ and P2), the effective critical 
stress is found to vary linearly with the mixing param- 
eter x. We have pointed out the origin of the apparent 
linear behaviour and argued that in general a non-linear 
variation is expected. 
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